Y = zeros(1,5000);
y0 = 0;
delt_T = 0.001;
T = 1;
uk_1 = 0;
uk = 0;
K=10;
u = 0;
for i=1:1: 5000
    if i==1000
        u=1;
    end
    uk_1 = uk;
    uk = u;
    yk = K* uk +K*T*((uk-uk_1)/delt_T);
    Y(i)=yk;
    U(i) =u;
end
t = 0: 0.001 : 4.999;
hold on
plot(t,Y);
plot(t,U);

